V 


■ 


Contract  Nonr  - 493(CX))»  fJR  010  104 
Principal  Investigator*  Robcr L G*  Parr 


Technical  Report  Number  5 
entitled 

ELECTRCNIC  STATES  OF  DlATO.tCC  MOLECULESi 
THE  OXYGEN  /.OLECULE 

by 

Fausto  G.  Fumi  and  Robert  G.  Parr 


CARNEGIE  INSTITUTE  OF  TECHT^LOGY 

f 

Department  of  Chemistry 


January  15»  19b3 


Pittsburgh  13»  Pennsylvania 


/ 
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A scheme  is  proposed  for  the  computation  of  electronic  energy 
levels  of  diatomic  molecules  with  fair  accuracy  and  minimum  labor*  The 
theoretical  ICAO  MO  method  including  configuration  Interaction  is  employ- 
ed. However*  the  number  of  configurations  considered  is  kept  small*  the 
energies  of  the  asymptotic  dissociation  products  for  the  various  states 
are  taken  from  atomic  data  in  the  manner  of  Parlser*  the  calculation  of 
the  interaction  energies  is  simplified  by  the  neglect  of  overlap  in  all 
terms  and  the  neglect  of  differential  overlap  in  the  electronic  repulsloh 
terms*  and  certain  core  parameters  are  fitted  empirically*  The  first 
application  is  made  to  the  oxygen  molecule*  using  one  empirical  para- 
meter* The  vertical  excitation  energies  from  the  ground  state  to  the 

and  ^2’  states  are  computed  as  functions  of  distance 
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between  1.16A  and  1.63A*  and  they  agree  with  the  observed  values  within 
0.2  ev;  excitation  energies  to  the  unobserved  ^2~*  ^Ay*^Ay  and  states 
are  also  computed*  The  results  are  in  substantial  agreement  with  a 
previous  more  involved  calculation  by  fAoffitt*  A justification  of  the 
proposed  scheme  is  presented  which  makes  use  of  the  orthonormallzed 
atomic  orbitals  of  Lowdin  * 
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1.-  The  availability  of  a large  amount  of  accurate  experimental  data 
on  the  electronic  levels  of  diatomic  molecules  provides  an  ideal  back- 
ground for  theoretical  work»  and  yet  few  quantitative  calculations  have 
been  made  on  diatomic  molecules  other  than  hydrogen  since  the  early 
applications  of  variational  techniques.  The  present  work  is  not  of 
comparable  rigor  but  rather  Involves  the  development  and  application 
of  a scheme  +o  compute  with  minimum  labor  and  fair  accuracy  electronic 
excitation  energies  in  diatomic  molecules  Involving  several  valence 
electrons*  without  necessarily  obtaining  a good  approximation  to  the 
ground  state . 

The  scheme  is  essentially  the  one  proposed  by  Pariser  and  Parr^ 
for  the  computation  of  the  energy  levels  of  complex  unsaturated  molecules. 
It  falls  within  the  framework  of  one* electron  theories  and  uses  the 
formalism  of  molecular  orbital  theory  in  the  approximation  in  which 
molecular  orbitals  (.';D)  are  expressed  as  linear  combinations  of  atomic 
orbitals  (ICAO)  . However*  we  consider  the  possibility  that  molecular 
states  might  not  be  properly  described  in  simple  configurational 
language*  where  a configuration  is  defined  by  assigning  electrons  to 
s — we  allow  for  cohfiguration  interaction.  The  choice  of  molecular 
orbital  language  is  suggested  by  the  availability  of  much  outstanding 
qualitative  work  by  f^lliken  and  others  on  the  interpretation  of  the 
electronic  spectra  of  diatomic  molecules  in  terms  of  VD  configurations. 

The  important  features  of  the  scheme^  are  the  obtalnment  of  the 
energies  of  the  asymptotic  dissociation  products  of  molecular  states 
from  atomic  data2»3,  a considerable  simplification  of  molecular  algebra^ 
effected  by  the  neglect  of  overlap  in  the  core  terms  and  the  neglect  of 
differential  overlap  in  the  electronic  repulsion  termsi  and  the  treatment 
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of  certain  molecular  core  parameters  as  empirical  quantities. 

For  a first  application  of  the  scheme  we  have  chosen  the  oxygen 
molecule.  Here  the  q*jantitatlve  treatment  is  sitrplified  by  the 
possibility  of  accounting  qualitatively  for  the  low  excited  states  on 
the  ir-electron  approximation,  /.toffitt^  and  .’teckler^  have  also  recently 
treated  02*  Mecklcr  was  interested  in  the  magnetic  properties  of  the 
ground  state  of  the  molecule  over  a wide  range  of  distance  and  used 
Gaussian  atomic  functions;  consequently  he  had  t.o  Include  a large 
amount  of  configuration  Interaction.  Moffitt  set  out  to  test  con- 
ventional IGAO  theory  with  configuration  Interaction*  and  he  arrived 
at  a preliminary  formulation  and  made  a first  application  of  his 
important  "atoms  in  molecules"  idea.  In  his  formulation  one  constructs 
molecular  wave  functions  of  proper  behaviour  at  infinite  distance;  this 
Is  not  required  in  the  present  formulation. 

2.-  The  qualitative  interpretation  of  the  low  excited  states  of  0^ 
involves  only  two  configurations: 

A;  [(KK)  ( <Tg2s)^(  CTu2s)2(  cTg2p)^](Try2p)^(7rg2p)2 

Bi  [(KK)  ( qTg2s)2{  C^2s)^(  q^2p)2j(Tr^2p)^(TTg2p)^ 

Indeed  one  can  show  by  using  the  vector  model  and  Pauli  principle  that 
the  A configuration  gives  the  states  '’2"*  2^  and  Ag»  while  the  B con- 
figuration gives  the  states  ^2y,  and  ^Ay.  These 

include  the  observed  low  excited  states  of  C>2* 

It  seems  reasonable  to  try  and  develop  a quantitative  treatment 
of  these  states  on  the  7r-electron  approximation^.  There  is  then  only 
one  configuration  to  consider  besides  A and  B: 

Cl  (^(KK)  { q-g2s)^{  <Ty2s)2(  <T^2p)2](Try2r  )2(Fg2p)'^ 
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The  choice  of  values  for  the  one^center  Integrals  (aa|aa)  and 
(aa|8sa$)  is  a crucial  point  in  the  scheme,  ftoffitt^  and  Pariser^  have 
independently  pointed  out  that  a major  source  of  error  in  ICAO  tJQ  cal- 
culations of  excitation  energies  using  Slater  AO' s is  the  inaccuracy 
of  the  energies  of  dissociation  products*  and  they  have  given  ways  of 
eliminating  this  error  which  involve  taking  data  from  atomic  spectroscopy 
and  using  atomic  wave  functions  only  to  compute  the  properly  molecular 
part  of  the  energy.  Y.'e  take  this  attitude  and  extend  the  method  of 
Parlser^* determining  the  values  of  (aa|aa)  and  (aa|asas)  from  an 
analysis  of  the  valence  states  of  O,  0"*“  and  0“  which  are  relevant  to  the 
molecular  states  being  considered.  The  identification  of  these  valence 
states  is  schematically  described  in  Table  VI;  lable  VII  gives  their 
location  in  energy  as  derived  from  a f\illiken-type  treatment^.  Values 
of  (aa|aa)  and  (aajasas)  are  obtained  by  a least  squares  method  from 
the  following  equations  relating  theoretical  and  experimental  energy 
changes  for  dissociation  processes: 

2 0(V2)  “ + 0”(Vx)f  (aa|aa)  + (aasl  aag)  «=  15.4 

2 0(V2)  “ 0+(Vi)  + O-(Vi),  (aa|aa)-(aaslea3)  « 13.7  (6) 

2 0(V2)  - O^C^Vj)  + O'(V^),  (aalaa)-3(aas|aas)»*^11.6 

For  the  two  center  Integrals  (aa|bb)  and  (aa|bgbs)  we  take  the 
values  given  by  Kopineck^®  for  Slater  2pn  AO's  of  effective  charge  4.55. 
"Ball"  values  for  these  integrals^  were  discarded  because  they  do  not 
satisfy  the  identity, 

(aajbb)  - (aalbgbs)  - 2(aas|bbsi 
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(7) 


- 7 - 

which  relates  integrals  over  Slater  orbitals*  and  this  identity  is  used 
in  obtaining  the  one-center  integrals  from  valence-state  data.  We  have 
also  found  that  lowering  the  values  of  (aalbb)  and  (aatbsbg)  as  suggested 
by  Parlser  and  Parr^  lessens  the  agreement  with  experiment  of  the  final 
excitation  energies. 

The  value  of  the  core  parameter  « is  not  needed  for  the  calculation 
of  the  excitation  energies.  We  choose  to  consider  p an  empirical 
quantity*  fitting  it  to  the  observed  excitation  energy  ^2“ 
obtained  from  ftorse  curves  for  these  two  states.  The  resulting  values 
of  p are  included  in  Table  V. 

5.-  Table  VIII  gives  the  electronic  excitation  energies  from  the 
ground  state  to  the  low  excited  states  of  O2  calculated  using  the  energy 
formulas  of  Tables  III  and  IV  and  the  values  for  the  integrals  given  in 

Table  V.  The  experimental  values  for  the  excitation  energies  quoted  in 

Table  VIII  are  obtained  by  plotting  Morse  curves  for  the  ground  state 
and  for  each  of  the  excited  states  using  the  experimental  values  for  the 
equilibrium  distance  re*  the  dissociation  energy  De  and  the  harmonic 
vibration  frequency  «Og7*ll,  the  states  and  the  resulting 

experimental  excitation  energies  are  reliable  to  a few  thousandths  of 
an  ev*  but  for  the  states  and  the  {precision  drops  to  a few 

hundredths  and  a few  tenths  of  an  ev*  respectively. 

Our  theoretical  results  agree  with  the  experimental  values  within 
0.2  ev.  For  the  state  ^JJ*  however*  the  agreement  is  somewhat  uncertain 
since  the  experimental  values  are  not  known  accurately;  anyhow*  it  is 
clear  that  for  this  state  the  ir-electron  approximation  breaks  down  at 
the  larger  distances  considered,  as  ;*>ffitt^  has  also  noted. 
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A conparison  of  our  calculated  values  with  those  of  ftoffitt 
reveals  very  good  general  agreement . The  approximation  of  zero 
differential  overlap  may  be  said  to  be  thus  vindicated . 

The  predicted  excitation  energies  to  the  non-observed  states  ^2^, 
and  should  be  somewhat  qualified.  We  regard  with  some 
confidence  the  predictions  for  the  states  ^2“  and  which  should  be 
stable  and  have  a value  for  of  the  order  of  1 ev.  The  predictions 
for  the  more  excited  states  ^Ay  and  ^2^,  on  the  other  hand*  probably 
have  only  a qualitative  value*  giving  the  general  location  and  the 
order  in  energy  of  these  two  states;  whether  they  should  be  stable 
remains  an  open  question. 

APPENDIX 

A justification  of  the  scheme 

The  features  of  the  scheme  that  need  to  be  justified  are  the 
neglect  of  overlap  in  all  energy  terms  and  the  neglect  of  differential 
overlap  In  the  electronic  repulsion  terms*  together  with  the  use  of 
Slater  AO* s to  compute  the  integrals  which  ere  not  assumed  equal  to  zero. 

Lowdin^^  has  proved  elegantly  that  his  orthonormalized  atomic 
orbitals  ("Lowdin  orbitals" ) ^2,  which  are  defined  to  satisfy  exactly 
the  condition  of  zero  overlap*  satisfy  approximately  the  condition  of 
zero  differential  overlap.  Here  we  shall  consider  in  detail  the  case 
of  diatomic  molecules*  deriving  the  relations  between  integrals  over 
Lowdln>  orbitals  and  integrals  over  Slater  orbitals . We  shall  see  that 
the  corrections  are  important  only  for  the  core  integrals  while  for  the 
electronic  repulsion  integrals  the  corrections  may  ordinarily  be 
neglected . 
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Consider  first  a homonuclear  diatomic  molecule  with  equivalent 
Slater  AO*s  Xs  and  Xb  atoms.  Denoting  by  and  the 

H 

corresponding  Lowdin  orbitals  we  have 


*b  • + AXb  • 


e^ere 

A ■ + 

B - - 
in  which  S 


2 - 2S2 


■l+^§^+  ••••  » 

1 

■»-2S+****» 

Then*  denoting  integrals  over  Lowdin  orbitals 


by  superscripts  X* 


» 0,  (9) 

* f^|§‘  » (10) 

/*‘  ■ ^ -/>-  , {ID 

Where  «<  ■ I 3^«dv  and  6»  are  the  Coulomb  and 

resonance  core  integrals.  The  formula  for  is  indeed  the  one  that 
Parlser  and  Parrl  propose  to  apply  in  computing  the  core  resonance 
intv'^gral  to  be  used  in  the  energy  formulas  whenever  one  does  not  want 
to  consider  this  integral  an  empirical  quantity.  If  we  use  the 
Mulliken®  approximate  formulas  for  the  hybrid  Coulomb-exchange  and  the 
exchange  integral  between  Xa  X b (taptaraount  to  considering  the 
first  terra  in  an  infiflte  series  expansion  in  powers  of  S^^)* 
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(aalab)  ■ ^f(aalaa)  (aa|bb}l* 

(12) 

(ab|ab)  ■ S (aa|ab)f 
we  obtain  also 

(aa\aa)^  ■ (aa|aa)'f [(aa|aa)-(aa|bb)]  (aa|aa)> 

(aa|bb)^  ■ (aa|bb)-(^^|jj)  [(aa|aa)-(aa|bb)J  ^ (aa|bb),  (13) 
(aa|ab)^  " (abjau)^  ■ O • 

Supprassion  of  the  Mulliken  approximation  would  cause  terms  in  to 
appear  in  these  last  formulas* 

If  each  atom  of  the  diatomic  molecule  carriesttwo  orthogonal 
Slater  A0*s  say  7C3  and  on  the  first  atom  and  and  ^ bs  on  the 

second*  as  in  02»  the  argument  can  be  generalized  without  difficulty. 
In  particular*  one  finds  that  the  integrals  over  iJ&wdin  orbitals 
satisfy  the  identity*  Eq*  (7)*  which  relates  Slater  orbitals* 
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TABLE  II.  Energy  formulas  in  terms  of  int.3^i.ils 
over  molecular  orbitals.®*^ 


I(^u) 

lUg) 

T / t * 

(irg^TTgJ 

J (TT^fTTg) 

%g(A) 

4 

2 

6 

1 

a 

325(C) 

2 

4 

1 

6 

8 

HaC 

125(A) 

4 

2 

6 

1 

8 

l2S(C) 

2 

4 

1 

6 

8 

”ac 

^Ag(A) 

4 

6 

1 

8 

2 

4 

1 

6 

8 

“ac 

3 

3 

3 

3 

9 

^2G(B) 

3 

3 

3 

3 

9 

^Au(B) 

3 

3 

3 

3 

9 

l2$(B) 

3 

3 

3 

3 

9 

^2G(B) 

3 

3 

3 

3 

9 

^Au(B) 

3 

3 

3 

3 

9 

* Except  for  the  Hac»  fhe  formulas  given  are  for  the  energies  of 
states  of  the  indicated  svwmetrles  for  the  indicated  configurations. 
The  Hyve  *re  interaction  elements  between  the  A and  C configurations. 


Definitions  of  integrals  are  l(^)=J^Hcnre^Pv,  J^i»04)“j^(l)0j(2) 
p/rig)0i{l)C^i(2)dv,  K(0i’0i)*50:  (l|^(2HeVri2)25vl%  (2)^v  and 
fi0i  0i0i)  - )0i(l)0i(2)to%i2)0i(l)0i(2)dv. 


TABLE  III*  Energy  forrsulas  in  terms  of  Integrals  over 

atomic  orbitals  — A and  C configurations.**^ 


i(*a|aa) 

^(aaiagas) 

J(aa| bb) 

yCaalbebj) 

Core  Contribution 

®*5(A) 

3 

23 

15 

19 

6^  * 2fL 

**5(c) 

3 

23 

15 

19 

6«  - 2p 

Hac 

-1 

3 

1 

-3 

0 

1*5  (A) 

7 

19 

19 

15 

+ 2p 

7 

19 

19 

15 

(f<  - 2f 

Hac 

3 

-1 

-3 

1 

0 

^Ag(A) 

5 

21 

17 

17 

6^  + 2p 

5 

21 

17 

17 

6»<  - 2(1 

Hac 

1 

1 

-1 

-1 

0 

* See  Table  II»  note  a* 

b Definitions  of  Integrals  are  (pqlrs)  • iXp(DXa(l)  (eVr,2)Xr(2)X«(2)  dv, 
■ jXaHcore^a<lv  and  i^bHcore3^a<^v,  where  Xp»Xq»X-  and  X 5 are*Slater 
AO*  8 from  the  setf  Xa,Xag»Xb,^bg  • 


TABLE  IV*  Energy  formulas  in  terms  of  integrals  over  atomic 
orbitals  --  B configuration.^'^ 


^(aajaa) 

^(aajasas) 

^(aajbb) 

|(aa|b,bg) 

Core  Contribution 

^u 

2 

10 

11 

7 

6^ 

u 

1 

- 13 

8 

8 

2 

10 

9 

9 

6^ 

5 

9 

8 

8 

ef* 

2 

10 

7 

11 

3 

11 

8 

8 

6«t 

^ See  Table  II>  note  a. 

^ See  Table  III*  note  b. 


TABLE  V.  Values  of  integrals  (ev) 


0 

Internuclear  Distance  (A) 

1 .163  1 .279  1 .396 

1. 

1.512 

1.628 

(aa|bb) 

14  .52 

11.213 

10 .350 

9.601 

8.947 

8.372 

(aalb^bg) 

12.62 

11.019 

10.216 

9.507 

8.880 

8.324 

5.381 

3.283 

2.390 

1.956 

1.689 

^ See  Table  III»  note  b*  for  definitions  of  integrals. 


TABLE  VI.  Valence  states  of  Ot  0“*^  and  0“  • 


State 

Schematic 

Representation* 

TT“ 

Energy  Formula^ 

> 

0k 

1 

o 

0 

0 

0 O 

^Jrnr'^Tnr 

4‘<(0') 

CM 

> 

0k 

o 

o o 

0 

o 
o o 

^7nr“  ■*■ 

3^(0) 

O'*’,  Vi 

{“ 

0 

0 o] 

Jw  + 

2«<(0'*') 

0'*‘,3vi 

0 

o 

^inr 

2«<(0'*‘) 

O’*-,  ^Vi 

c 

o 

c “) 

‘^Tnr'*'  ^mr  ■*■ 

* An  o in  column  under  f j indicates  occupancy  of  f 

by  an  electron  of  spin  . Brackets  denote  alternativCf  parenthese 
coupled  situations. 

^ *=<(0)  and  ®<{0^)  denote  the  core  energies  of  a 2pir  electron  in 

C“»  O and  0"^  respectively.  is  the  Coulomb  integral*  Kyrn  "the 

exchange  Integral  between  two  2pw  electrons  in  0~*  0 or  0'*’»  as  the 
case  may  be. 


rm 


TABLE  VII.  Energies  of  spectroscopic  and  valence  states  of  0> 


O’*"  and  0“  (ev.)  . 


State 

Energy 

Details  of  Energy  Calculation®'^ 

0”^*  ^/l»0”*’*  ^ 

18.6 

18.6  * 13.6't-5.0*  5.0  from  spectral  data  for  O'** 

0+* 

16.9 

16.9  * 13.6t3.3»  3.3  from  spectral  data  for  O'*" 

0*,h/i 

14.7 

14.7  « 13.6+1.1*  1.1  « ^^3^ 

0'*‘*^S 

13 .6 

ionization  potential  of  0 

> 

O 

0.5 

0.5  • 4 ^^2^^  * spectral  data  for  0 

0 *3p 

0.0 

arbitrary  energy  zero 

0“*Vi*0-*2p 

-2.2 

negative  electron  affinity  of  C 

* Spectral  data  for  0 and  0*  and  the  ionization  potential  and  electron 
affinity  for  O are  taken  from  Landolt-Bornstein  Tables*  Vol.  I*  1950. 

^ In  the  two  cases  in  which  averages  of  spectroscopic  states  are  taken* 

I the  averages  are  formed  in  such  a way  as  to  give  energy  values  in  agree- 
I ment  with  those  quoted  by  f.k>ffitt  (reference  5)  . Alternative  averaging 

I processes  are  possible  which  give  slightly  different  results. 

I 

( 


♦ 

f 

} 


TABLE  VIII*  Elttctronic  •xcltatlon  •nergles  in  O2  (ev) *> 


Int«rnucl«ax  Distance  (A)^ 


1.163 

1.279 

1.396 

1.512 

1.628 

% 

calculated 

0.96 

0.86 

0.72 

0.59 

0.48 

observed 

1.01 

0.94 

0.84 

0.73 

0.62 

calculated 

1.90 

1.61 

1.30 

1.04 

0.83 

9 

observed 

1.70 

1.55 

1.37 

1.17 

0.97 

calculated 

observed 

calculated 

observed 

10.39 

5.85 

3.80 

2.75 

2.09 

10 .59 

5.98 

3.90 

2.81 

2.14 

3j+ 

calculated 

10 .78 

6.12 

3.99 

2.88 

2.18 

u 

observed 

10 .78 

6.12 

3.99 

2.88 

2.18 

c 1 

calculated 

11.14 

7.37 

6.01 

5.57 

5.46 

observed 

11.26 

7.94 

5.91 

4.65 

3.83 

\ 

calculated 

observed 

calculated 

observed 

13.16 

9.84 

7.81 

6.55 

5.73 

u 

15 .06 

11.74 

9.71 

8.45 

7.63 

* Excitation  energies  are  relative  to  the  ground  state. 

^ The  equilibrium  intemuclear  distance  for  the  ground  state 
is  1.208A. 


